Conceptual Problems in Land Surface Data Assimilation

A land data assimilation system (LDAS) merges observations (or satellite retrievals) of land surface hydrological conditions, including soil moisture, snow, and terrestrial water storage (TWS), into a numerical model of land surface processes. In theory, the output from such a system is superior to estimates based on the observations or the model alone, thereby enhancing our ability to understand, monitor, and predict key elements of the terrestrial water cycle. In practice, however, several conceptual problems can interfere with realizing the potential improvements from data assimilation. Of particular concern is the frequent mismatch between the assimilated observations and the land surface model variables of interest. The seminar will discuss recent research with the ensemble-based NASA GEOS-S LDAS to address various aspects of this mismatch. These aspects include (i) the assimilation of coarse-scale observations into higher-resolution land surface models, (ii) the partitioning of satellite observations (such as TWS retrievals) into their constituent water cycle components, (iii) the forward modeling of microwave brightness temperatures over land for radiance-based land surface data aSSimilation, and (iv) the selection of the most relevant types of observations for the analysis of a specific water cycle variable (such as root zone soil moisture). At its core, the solution to the above challenges involves the careful construction of an observation operator that maps from the land surface model variables of interest to the space of the assimilated observations

Observation Error Specification for Near-Surface Soil Moisture Assimilation: Does It Matter?

The observation error variances assumed in the data assimilation (DA) of remotely sensed near-surface soil moisture observations are known to have unrealistic spatial patterns. The temporal variability in the errors is also typically neglected. Soil moisture errors cannot be directly calculated at large scales, and the globally distributed errors are unknown. We apply several approaches to estimating the observation errors (R matrix) for near-surface soil moisture DA experiments, to establish whether the assumed errors are important to the DA performance, and hence whether uncertainty in the observation error specification is a likely source of uncertainty in soil moisture DA output

Observation-Corrected Precipitation Estimates in GEOS-5

Several GEOS-5 applications, including the GEOS-5 seasonal forecasting system and the MERRA-Land data product, rely on global precipitation data that have been corrected with satellite and or gauge-based precipitation observations. This document describes the methodology used to generate the corrected precipitation estimates and their use in GEOS-5 applications. The corrected precipitation estimates are derived by disaggregating publicly available, observationally based, global precipitation products from daily or pentad totals to hourly accumulations using background precipitation estimates from the GEOS-5 atmospheric data assimilation system. Depending on the specific combination of the observational precipitation product and the GEOS-5 background estimates, the observational product may also be downscaled in space. The resulting corrected precipitation data product is at the finer temporal and spatial resolution of the GEOS-5 background and matches the observed precipitation at the coarser scale of the observational product, separately for each day (or pentad) and each grid cell

Assessment of MERRA-2 Land Surface Energy Flux Estimates

In MERRA-2, observed precipitation is inserted in place of model-generated precipitation at the land surface. The use of observed precipitation was originally developed for MERRA-Land(a land-only replay of MERRA with model-generated precipitation replaced with observations).Previously shown that the land hydrology in MERRA-2 and MERRA-Land is better than MERRA. We test whether the improved land surface hydrology in MERRA-2 leads to the expected improvements in the land surface energy fluxes and 2 m air temperatures (T2m)

Spring Hydrology Determines Summer Net Carbon Uptake in Northern Ecosystems

Increased photosynthetic activity and enhanced seasonal CO2 exchange of northern ecosystems have been observed from a variety of sources including satellite vegetation indices (such as the Normalized Difference Vegetation Index; NDVI) and atmospheric CO2 measurements. Most of these changes have been attributed to strong warming trends in the northern high latitudes (greater than or equal to 50N). Here we analyze the interannual variation of summer net carbon uptake derived from atmospheric CO2 measurements and satellite NDVI in relation to surface meteorology from regional observational records. We find that increases in spring precipitation and snow pack promote summer net carbon uptake of northern ecosystems independent of air temperature effects. However, satellite NDVI measurements still show an overall benefit of summer photosynthetic activity from regional warming and limited impact of spring precipitation. This discrepancy is attributed to a similar response of photosynthesis and respiration to warming and thus reduced sensitivity of net ecosystem carbon uptake to temperature. Further analysis of boreal tower eddy covariance CO2 flux measurements indicates that summer net carbon uptake is positively correlated with early growing-season surface soil moisture, which is also strongly affected by spring precipitation and snow pack based on analysis of satellite soil moisture retrievals. This is attributed to strong regulation of spring hydrology on soil respiration in relatively wet boreal and arctic ecosystems. These results document the important role of spring hydrology in determining summer net carbon uptake and contrast with prevailing assumptions of dominant cold temperature limitations to high-latitude ecosystems. Our results indicate potentially stronger coupling of boreal/arctic water and carbon cycles with continued regional warming trends

Inferring Land Surface Model Parameters for the Assimilation of Satellite-Based L-Band Brightness Temperature Observations into a Soil Moisture Analysis System

The Soil Moisture and Ocean Salinity (SMOS) satellite mission provides global measurements of L-band brightness temperatures at horizontal and vertical polarization and a variety of incidence angles that are sensitive to moisture and temperature conditions in the top few centimeters of the soil. These L-band observations can therefore be assimilated into a land surface model to obtain surface and root zone soil moisture estimates. As part of the observation operator, such an assimilation system requires a radiative transfer model (RTM) that converts geophysical fields (including soil moisture and soil temperature) into modeled L-band brightness temperatures. At the global scale, the RTM parameters and the climatological soil moisture conditions are still poorly known. Using look-up tables from the literature to estimate the RTM parameters usually results in modeled L-band brightness temperatures that are strongly biased against the SMOS observations, with biases varying regionally and seasonally. Such biases must be addressed within the land data assimilation system. In this presentation, the estimation of the RTM parameters is discussed for the NASA GEOS-5 land data assimilation system, which is based on the ensemble Kalman filter (EnKF) and the Catchment land surface model. In the GEOS-5 land data assimilation system, soil moisture and brightness temperature biases are addressed in three stages. First, the global soil properties and soil hydraulic parameters that are used in the Catchment model were revised to minimize the bias in the modeled soil moisture, as verified against available in situ soil moisture measurements. Second, key parameters of the "tau-omega" RTM were calibrated prior to data assimilation using an objective function that minimizes the climatological differences between the modeled L-band brightness temperatures and the corresponding SMOS observations. Calibrated parameters include soil roughness parameters, vegetation structure parameters, and the single scattering albedo. After this climatological calibration, the modeling system can provide L-band brightness temperatures with a global mean absolute bias of less than 10K against SMOS observations, across multiple incidence angles and for horizontal and vertical polarization. Third, seasonal and regional variations in the residual biases are addressed by estimating the vegetation optical depth through state augmentation during the assimilation of the L-band brightness temperatures. This strategy, tested here with SMOS data, is part of the baseline approach for the Level 4 Surface and Root Zone Soil Moisture data product from the planned Soil Moisture Active Passive (SMAP) satellite mission

Contribution of Soil Moisture Information to Streamflow Prediction in the Snowmelt Season: A Continental-Scale Analysis

In areas dominated by winter snowcover, the prediction of streamflow during the snowmelt season may benefit from three pieces of information: (i) the accurate prediction of weather variability (precipitation, etc.) leading up to and during the snowmelt season, (ii) estimates of the amount of snow present during the winter season, and (iii) estimates of the amount of soil moisture underlying the snowpack during the winter season. The importance of accurate meteorological predictions and wintertime snow estimates is obvious. The contribution of soil moisture to streamflow prediction is more subtle yet potentially very important. If the soil is dry below the snowpack, a significant fraction of the snowmelt may be lost to streamflow and potential reservoir storage, since it may infiltrate the soil instead for later evaporation. Such evaporative losses are presumably smaller if the soil below the snowpack is wet. In this paper, we use a state-of-the-art land surface model to quantify the contribution of wintertime snow and soil moisture information -- both together and separately -- to skill in forecasting springtime streamflow. We find that soil moisture information indeed contributes significantly to streamflow prediction skill

Variational assimilation of remote sensing data for land surface hydrologic applications

Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Civil and Environmental Engineering, 2000.Includes bibliographical references (p. 283-192).Soil moisture plays a major role in the global hydrologic cycle. Most importantly, soil moisture controls the partitioning of available energy at the land surface into latent and sensible heat fluxes. We investigate the feasibility of estimating large-scale soil moisture profiles and related land surface variables from low-frequency (L-band) passive microwave remote sensing observations using weak-constraint variational data assimilation. We extend the iterated indirect representer method, which is based on the adjoint of the hydrologic model, to suit our application. The four-dimensional (space and time) data assimilation algorithm takes into account model and measurement uncertainties and provides optimal estimates by implicitly propagating the full error covariances. Explicit expressions for the posterior error covariances are also derived. We achieve a dynamically consistent interpolation and extrapolation of the remote sensing data in space and time, or equivalently, a continuous update of the model predictions from the data. Our hydrologic model of water and energy exchange at the land surface is expressly designed for data assimilation. It captures the key physical processes while remaining computationally efficient. The assimilation algorithm is tested with a series of experiments using synthetically generated system and measurement noise. In a realistic environment based on the Southern Great Plains 1997 (SGP97) hydrology experiment, we assess the performance of the algorithm under ideal and non ideal assimilation conditions. Specifically, we address five topics which are crucial to the design of an operational soil moisture assimilation system. (1) We show that soil moisture can be satisfactorily estimated at scales finer than the resolution of the brightness images (downscaling), provided sufficiently accurate fine-scale model inputs are available. (2) The satellite repeat cycle should be shorter than the average interstorm period. (3) The loss of optimality by using shorter assimilation intervals is offset by a substantial gain in computational efficiency. (4) Soil moisture can be satisfactorily estimated even if quantitative precipitation data are not available. (5) The assimilation algorithm is only weakly sensitive to inaccurate specification of the soil hydraulic properties. In summary, we demonstrate the feasibility of large-scale land surface data assimilation from passive microwave observations.by Rolf H. Reichle.Ph.D

Spring hydrology determines summer net carbon uptake in northern ecosystems

Increased photosynthetic activity and enhanced seasonal CO2 exchange of northern ecosystems have been observed from a variety of sources including satellite vegetation indices (such as the normalized difference vegetation index; NDVI) and atmospheric CO2 measurements. Most of these changes have been attributed to strong warming trends in the northern high latitudes (50° N). Here we analyze the interannual variation of summer net carbon uptake derived from atmospheric CO2 measurements and satellite NDVI in relation to surface meteorology from regional observational records. We find that increases in spring precipitation and snow pack promote summer net carbon uptake of northern ecosystems independent of air temperature effects. However, satellite NDVI measurements still show an overall benefit of summer photosynthetic activity from regional warming and limited impact of spring precipitation. This discrepancy is attributed to a similar response of photosynthesis and respiration to warming and thus reduced sensitivity of net ecosystem carbon uptake to temperature. Further analysis of boreal tower eddy covariance CO2 flux measurements indicates that summer net carbon uptake is positively correlated with early growing-season surface soil moisture, which is also strongly affected by spring precipitation and snow pack based on analysis of satellite soil moisture retrievals. This is attributed to strong regulation of spring hydrology on soil respiration in relatively wet boreal and arctic ecosystems. These results document the important role of spring hydrology in determining summer net carbon uptake and contrast with prevailing assumptions of dominant cold temperature limitations to high-latitude ecosystems. Our results indicate potentially stronger coupling of boreal/arctic water and carbon cycles with continued regional warming trends

Multiscale assimilation of Advanced Microwave Scanning Radiometer-EOS snow water equivalent and Moderate Resolution Imaging Spectroradiometer snow cover fraction observations in northern Colorado

Eight years (2002–2010) of Advanced Microwave Scanning Radiometer–EOS (AMSR-E) snow water equivalent (SWE) retrievals and Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) observations are assimilated separately or jointly into the Noah land surface model over a domain in Northern Colorado. A multiscale ensemble Kalman filter (EnKF) is used, supplemented with a rule-based update. The satellite data are either left unscaled or are scaled for anomaly assimilation. The results are validated against in situ observations at 14 high-elevation Snowpack Telemetry (SNOTEL) sites with typically deep snow and at 4 lower-elevation Cooperative Observer Program (COOP) sites. Assimilation of coarse-scale AMSR-E SWE and fine-scale MODIS SCF observations both result in realistic spatial SWE patterns. At COOP sites with shallow snowpacks, AMSR-E SWE and MODIS SCF data assimilation are beneficial separately, and joint SWE and SCF assimilation yields significantly improved root-mean-square error and correlation values for scaled and unscaled data assimilation. In areas of deep snow where the SNOTEL sites are located, however, AMSR-E retrievals are typically biased low and assimilation without prior scaling leads to degraded SWE estimates. Anomaly SWE assimilation could not improve the interannual SWE variations in the assimilation results because the AMSR-E retrievals lack realistic interannual variability in deep snowpacks. SCF assimilation has only a marginal impact at the SNOTEL locations because these sites experience extended periods of near-complete snow cover. Across all sites, SCF assimilation improves the timing of the onset of the snow season but without a net improvement of SWE amounts